Intrinsic Spin Hall Effect in s-wave Superconducting State: 
Analysis of Rashba Model 
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A general expression for the spin Hall conductivity (SHC) in the s-wave superconducting state 
at finite temperatures is derived. Based on the expression, we study the SHC in a two-dimensional 
electron gas model in the presence of Rashba spin-orbit interaction (SOI). The SHC is zero in the 
normal state, whereas it takes a large negative value as soon as the superconductivity occurs, due to 
the change in the quasiparticle contributions. Since this remarkable behavior is independent of the 
strength of the SOI, it will be widely observed in thin films of superconductors with surface-induced 
Rashba SOI, or in various non-centrosymmetric superconductors. 
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The spin Hall effect (SHE), which is a phenomenon 
that an electric field E induces a transverse spin cur- 
rent J s , has attracted considerable attention. The in- 
trinsic SHE, which is independent of impurity scatter- 
ing, was first predicted in semiconductors [1-4]. Later, 
it was revealed that the large intrinsic SHE emerges in 
transition metals and their compounds, since the com- 
plex ci-orbital wavefunction induced by spin-orbit inter- 
action (SOI) gives rise to the "orbital Berry phase" that 
works as a spin-dependent effective magnetic field [5-8]. 
In several metals, the spin Hall conductivity (SHC, ctsh) 
had been experimentally determined by observing the in- 
verse SHE signal [9-11]. The observed SHC in Pt exceeds 
200 He' 1 -Q- 1 cur 1 [10], and the SHCs in Nb and Mo are 
negative [12]. These results can be explained in terms of 
the intrinsic SHE [8, 13], suggesting the importance of 
the intrinsic mechanism in transition metals. 

In spintronics, superconductivity is widely used to con- 
trol the spin state. In this sense, a natural question is 
whether and how the SHE emerges in the superconduct- 
ing state. Although DC electric field E cannot exists in 
bulk superconductors, the SHC in superconductors is ob- 
servable as follows: For example, the Hall voltage due to 
inverse SHE should appear in the superconducting tun- 
neling junction, thanks to the charge imbalance effect 
[14]. Moreover, the gradient of temperature will induces 
the transverse spin current in a bulk superconductors 
when the SHC is finite: The observed spin Nernst con- 
ductivity «sh = Jy/(~^xT) [15] would be related to the 
SHC by the Mott relation ctsh oc Tcrg H (ep), as discussed 
in the intrinsic anomalous Hall effect [16]. In addition, 
the SHC in type-II superconductors will be measurable in 
the mixed state under the magnetic field, since the resis- 
tivity is nonzero [17]. Therein, SHC will not be affected 
by moving vortices since they do not convey spin. 

In this letter, we present a general expression for the 
intrinsic SHC in the superconducting state. This expres- 
sion indicates that the current vertex correction (CVC), 
which is a consequence of the conservation laws in the 
field theory, can cause a significant change in the SHC in 



the superconducting state. Based on the expression, we 
analyze the two-dimensional electron gas (2DEG) model 
with Rashba SOI as a typical spin Hall system, and re- 
veal that a giant SHC emerges in the superconducting 
state due to the CVC, independently of the strength of 
the SOI. This phenomenon will be observed in thin film 
superconductors with the aid of surface-induced Rashba 
SOI [18], or in various non-centrosymmetric supercon- 
ductors such as CePtaSi, L^PtaB, MgSi, and Y2C3. 

The Rashba 2DEG model describes the electrons in a 
semiconductor inversion layer. It had been studied inten- 
sively not only as the issue of semiconductor spintronics, 
but also as a model for non-centrosymmetric supercon- 
ductor [19]. The SHC in Rashba 2DEG model was first 
studied by Sinova et al [1]. They found that the SHC 
takes a universal value —e/8n in the ballistic regime. (— e 
is the electron charge.) Later, Inoue et al [3] had shown 
that the SHC vanishes identically in the diffusive regime 
due to the CVC induced by short-range impurities. This 

fact is also explained by the relation S x = i[H, S x ] cx 
in the Rashba model [4, 20], where is the spin cur- 
rent. However, this relation does not hold in the super- 
conducting state. For this reason, the SHC in Rashba 
2DEG model shows a considerably large value just below 
T c unless the SOI is zero, as we will show below. 

The s-wave BCS- Rashba 2DEG model is expressed by 
the following 4x4 form in the Nambu representation [19] : 



/i k — iAa y 



-h°* 



(1) 



where &i (i = x,y,z) is the Pauli matrix, and </> k = 
( c kt, cL , c_k|i c -ki); c kcr is an annihilation operator of 



electron. /i k = e kl + ^{& x k y — cr y k x ) is the 2 x 2 k- 
linear Rashba 2DEG Hamiltonian in the normal state, 
where A is the SOI parameter and ek = k 2 /2m — /i 
is the chemical potential. A is the superconducting gap; 
we assume that < A -C fi hereafter. The quasiparticle 
spectrum is given by Ek,± = \J (ek ± Afc) 2 + A 2 . We also 
we assume that the spin-orbit splitting, Aso = 2AA;f, is 
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much smaller than [i. In this case, the triplet pairing 
induced by the Rashba SOI is negligible [19]. 

Below, we study the SHC in the presence of short-range 
nonmagnetic impurities, by which T c is unchanged due to 
Anderson's theorem [21]. In the Nambu representation 
[21], a single impurity potential term is given by 



>k', 



(2) 



k,k' 



where I is the impurity potential. Here and hereafter, 
we define the following 4x4 unitary matrices: 



To = 
f 2 = 



1 

o -i 

a z 
a z 



T 3 = 



t x 
-a, 

a z 
a z 



(3) 



The Green function in the Nambu representation is 
Gk(oj) = (ojI — i?k) _1 [21], and the retarded and ad- 
vanced Green function are G R {io) = G^(uj + ij) and 
G A (u>) = G^(uj — In the Born approximation, the 
damping rate is given by 7 = n im p TrI 2 N(uj), where a ^ ig given as 



The second term, AA^ = n imp I 2 £ k f G R A^G A f , rep- 
resents the CVC for /-term. As shown in Ref. [3] , the fac- 
tor n imp in eq. (7) cancels with ^2^.G R G A ~ 0(7 _1 ) ~ 
^( n imp)- Therefore, the CVC for /-term is important 
even in the clean limit (rii mp <C 1). We will show this 
fact explicitly in later calculation. 

On the other hand, CVC for //-term is negligible in 
the clean limit: The charge CVC for //-term, AA^ 77 , 

is given by eq. (7) by replacing G A with G R . Since 
J2k G R G R is regular for 7 -> +0, AA£ 77 ~ 0(n imp ), 
which is negligible in the clean limit. Moreover, AA^ 77 is 
related to the self-energy E k through the Ward identity. 
Since £ k is k-independent in the SCBA for short-range 
impurities, AA^ 77 vanishes in the present study. Thus, 
bare currents enter into eq. (5). 

By following Ref. [8], we rewrite eqs. (4) and (5) in 
the band-diagonal basis. In this basis, (UkG R Uk)i,m = 
[G R ]i, m = Si, m /(u - El), where £ k = ±E K± and L> k 
is the unitary matrix for the change of basis. Then, we 
can perform (^-integrations using the relationships Im(z— 



= ttS(z) and lm(z— 



-ttS'(z). For example, 



N(uj) = (m/27r)Re{l/ v /l - (A/w) 2 } is the density of 
states (DOS) per spin. 

Since the SHC is independent of the Meissner current, 
it will be insensitive to the electric field frequency oj for 
\lo\ <C 7o- In the linear response theory [8, 22], the DC 
SHC is expressed in the 4x4 Nambu representation as 
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where [JJ] m .i = hii{[Jf ] TO ,i[J° ];, TO }. The first (second) 
term corresponds to Og 7 ^ (c|h) in Refs. [8]. As a result, 
°sh is given by the summation of Cg^ 770 an d <y R ^: 



du 



T I+IIa 
7 SH 



(5) 
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[JA A] m , t 

[JJ]m,l 

2jEf - Elf 



du 



(8) 



[/(£<)- /(££")], (9) 



where f(u) — (e w / T + Jy is the bare charge 

current operator, which is given by — Sh k+eA /6A y \a=o 
((5/iI k+eA /5A y |A=o) for particle (hole) channel [22]. 
jj 1 = { j£, s z }/(— 2e) is the bare spin current operator. 
In the present model, they are given by 



where [JAA] TO: ; = Im{[jf ] m>; [AA°]; iTO }. Now, we can 
calculate the SHC in general superconductors using Eqs. 
(8) and (9), which had not been derived previously. 



J, = ~e -*l + ATi 



in 



2m 



(6) 




where — e\T\ in jy is called the anomalous velocity [1, 3], 
which is essential for the SHE. 

A^ in eq. (4) is the total charge current dressed by 
the CVC. When the elastic scattering dominates the in- 
clastic scattering, the CVC is derived from the following 
Bethe-Salpeter equation in the self-consistent Born ap- 
proximation (SCBA): 



A c 

v 



j: 



AA; 



(7) 



-kr 



FIG. 1: An example of terms for <r| H that is finite under T c . 
This term is proportional to A tx V / T~— T/T c . 

a s~H IIa m ec l- (8) represents the quasiparticle contribu- 
tion at the Fermi level, and in eq. (9) represents the 
Berry curvature term that is rewritten as the summation 
of the Berry curvature of Bloch electrons inside the Fermi 
sea. In the normal state in bulk transition metals, CVC 
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due to local impurities almost vanishes in the clean limit 
since the main matrix elements of current operators are 
odd functions of k. Therefore, al^ IIa w and the total 
SHC is approximately given by <r|^ for 7 — > [8]. In 
the Rashba 2DEG model, in contrast, crg^ //a = e/8n in 
the normal state because of the CVC for /-term. Since 
ct|jj = — e/8ir, the total SHC vanishes identically in the 
normal state [3]. 

Interestingly, such a balance between the quasiparti- 
cle term and the Berry curvature term in the normal 
state is drastically changed when the superconductivity 
(A 7^ 0) sets in: OgH 77a vanishes at T = because of the 
factor —df/duj. Thus, Cgj^ //a should show a prominent 
change below T c in spin Hall systems where the CVC is 
significant: such examples are metals with Rashba- type 
or Dresselhaus-type SOI [3, 4] and doped graphene [23]. 
On the other hand, <t|jj is finite even at T — unless A is 
much larger than the SOI splitting. Therefore, the SHC 
can show a nontrivial temperature dependence below T c 
due to the imbalance between cr^ IIa and Ogjj. 

Hereafter, we calculate both CgH //a and in the 
BCS-Rashba 2DEG model, and show that the SHC takes 
a large value just below T c . According to eq. (7), the 
lowest order CVC for AA^ is given as 



for x — 0; A. c 



AAf V) = 



''imp 



k 



(10) 



where A{u) = T G H JyG A T . Hereafter, we assume the 
damping rate in the normal state, 70 = rtj mp / 2 m/2, is 
the smallest parameter. Considering that the pole of 
G R is (j = ±E± - i-f and the relation (z 2 + 7 2 ) -1 = 
(tt/j)8(z), the expression of A(u>) is greatly simplified 
for 70 -C 1. After performing the integration J v d4>k/^ 
(4>k = t&n~ 1 (k v /k x )), A(uj) is simply given as 



A(u 



-e ■ 7r 

87 



Q = ±l 



AU(|c 



E, 



^(fi-a^fs) 



where x — A/uj. As a result, eq. (10) becomes 

AA^V) = -e(\/2)(-f 1+ xf 3 ), (11) 

where the relation 7 = 7mi mp / 2 N(u) is considered. 
The second lowest order CVC is given by cq. (10) 

by replacing Jy with AA^ . After taking the k- 
summation, it is simply obtained as 

AAf 2 V) = (l-x 2 )/2-Akf\u). (12) 
Thus, the total CVC, AA£ = J^Zi AA £ ( °. is § iven as 



<Tg H vanishes in the normal state since the off-diagonal 
anomalous velocity cancels out in the total charge current 



AAC 



-e(k y /m)l [3]. However, 



this cancellation does not occur in the superconducting 
state, and therefore alp ^ below T c . 

Now, we derive Ogjj that represents the quasiparticlc 
contribution. In a general basis, eq. (8) is rewritten as 



(7, 



I+IL 
SH 



-?/£(-£) 



B(u) 



(14) 



where B{uj) = ±Tr ^J^G R AA^G A . Since B(u) is sim- 
plified as (e/32m(l + a; 2 )) E Q =±i( e k + a\k) 2 5{\uj\ - E a ) 
for 7 <C 1, we obtain the following expression: 



T I+IIa 
'SH 



X(A,T) 



8tt 



X(A,T), 



(15) 



j —00 \ 



df] 
dtu) 



Ve 2 +A 2 



e 2 + 2A 2 ' 



It should be stressed that (y I s ^ IIa is independent of A(^ 
0). The asymptotic behavior of A (A, T) for A > is 

X(A,T) w 1 -7=777 for A < T (T - T c ), (16) 




e~ A / T for A> T (T<T C ). (17) 



In the same way, we derive that Og 7 jf = (e/8n)Y(A,T), 
where Y = de(—df/duj)^ (i i +A i is the Yosida func- 
tion. Since Y = 1 - 0((A/T c ) 2 ) for T ~ T c , cr| H w 
(-e/16\/2)(A/T). Fi gure 1 shows an example of terms 
for Og H that is of order O(A) after k, w-integrations. It 
represents the spin current due to the triplet (k f, — k f) 
particle-particle excitation induced by the Rashba SOI 
and the impurity scattering. 

Next, we discuss Og^, which is called the "Berry cur- 
vature term" [8]. It is caused by electrons in the Fermi 
sea that satisfy E l k ■ E^ < 0. In the normal state, the 
k-summation in eq. (9) is restricted to /cp- < |k| < &f+, 
where k^± are two Fermi momenta [8]. This restric- 
tion approximately holds in the superconducting state 
in the present model. Moreover, [Jf ] m ,z = [</^]m,z = 
for (E^jEU = (E^ a , —E^a), since there is no number- 
nonconserving elements in eq. (6) [22]. For these reasons, 
in contrast to Og^ 77 ", Og^ is insensitive to T and A if 
(T, A) <C Ago, which will be verified later by performing 
numerical calculation of eq. (9). 

Now, we show numerical results for m = fi = 1, under 
the condition that A, Ago 3> 7o- Figure 2 (a) shows the 
A-dependence of the SHC for A = 0.1 and T = 0.01. The 
quasiparticle term crg^ 77 " and the Berry curvature term 
(Tgn are given by eqs. (15) and (9), respectively To verify 
the correctness of cq. (15), we have also calculated the 
CVC by solving the Bethe-Salpeter eq. (7) numerically 
for 70 = 0.001, and derived <7g£ 77a using eq. (14). The 
obtained result is shown by square dots. Since cr| H = 
-(e/16V2)(A/T c ) + 0(A 2 ) and a£ H = 0(A 2 ), a| H gives 
the increment in |<tsh| just below T c . 
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SOI shows a prominent increment in the superconduct- 
ing state. Similar drastic change in the SHC will be ob- 
served in cases of the fc 3 -type Rashba or Dresselhaus SOI, 
since the CVC for the anomalous velocity is large in these 
cases. The present study opens the way to distinguish 
between the quasiparticle contribution (eq. (8)) and the 
Berry curvature contribution (eq. (9)), which had been 
desired to understand the mechanism of intrinsic SHE. 

In summary, we have presented the first study of the in- 
trinsic SHE in the superconducting state. In the Rashba 
2DEG model, the SHC changes from zero to a large nega- 
tive value just below T c unless A — 0, due to the change in 
the quasiparticle contribution Cg^ 77 ". This phenomenon 
will be observed in thin film superconductors with the aid 
of surface-induced Rashba SOI, in non-centrosymmetric 
superconductors, and in semiconductors using the super- 
conducting proximity effect. Moreover, the SHC in su- 
perconductors will also be observed using the AC mea- 
surement, since a xx {u>) is finite for u> ^ whereas the 
intrinsic SHC is independent of u> for < ui < -f~ 1 [3]. 

We are grateful to Y.Otani, T. Kimura, and E. Saitoh 
for valuable discussions. 



FIG. 2: (color online) (a) A-dependence of the SHC; com- 
parison of I-, Ila-, and 776-terms. (b) t = T/T c -dependence 
of the total SHC, ctsh, for T c = 0.2 ~ 0.005. 

Figure 2 (b) shows the T-dependence of the SHC 
for A = 0.05, assuming that A = Ao^l - T/T c and 
Ao = 1.8T C in accord with the weak-coupling BCS the- 
ory. It is noteworthy that both <r| H and crg^ /7a are unique 
functions of t = T/T c , and independent of T c and Aso 
c7g H takes a minimum value o-^ ax ~ —0.48 at t as 0.78. 
In the case of A < A so (i.e., T c < A so ), <J§^ as -e/8ir 
for < T < T c since the modification of electronic states 
for k-p- < |k| < /«f+ is limited, as discussed above. For 
this reason, ct sh as -e/8ir ■ (X(A,T) - 1) for A « A so , 
as shown in Fig. 2 (b). ct S h(T = 0) = <r££(T = 0) de- 
creases with Ao since the formation of singlet paring pre- 
vents the spin current. In the opposite case, Ao 3> Aso, 
the relation ctsh ~ 0g H is recognized, which means that 
°sh = Csif + °sh 1S verv small. Thus, SHC shows a 
striking enhancement just below T c for any value of A. 

In the present study, we consider that the conduction 
electrons are delocalized for fc F _ < |k| < k F+ , and they 
contribute to Og^. In conventional insulators, in con- 
trast, <Tg^ is expected to vanish since the electrons near 
the band edge are localized [24, 25]. (Note that local- 
ization cannot be described in the SCBA.) Interestingly, 
recent theoretical and experimental efforts have revealed 
that CTg^ takes a finite value in "topological insulators" 
such as graphene [24] and HgTe [25], owing to the delo- 
calized nature of massless Dirac fermions. 

We have shows that the SHE driven by fc-linear Rashba 
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